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We report in-plane 63 Cu Nuclear Magnetic Resonance measurements for a series of fully enriched 
(Cao.iLao.9)(Bai.65Lao.35)Cu30j / powder samples, which belong to the YBCO family, but the doping 
could vary from very underdoped to extremely overdoped. From these measurements we determine 
the average nuclear quadrupole resonance frequency vq, and its second moment Avq, both set by 
the in-plane hole density n, as a function of oxygen level y. We find that in the overdoped side 
n is saturated, but Avq rapidly increases with increasing y. The relevance of these results to the 
increasing penetration depth in overdoped cuprates is discussed. 



A common feature in the phase diagram of all cuprate 
high temperature superconductors (HTSC) is the dome- 
shaped curve of the T c versus doping, divided by an 
optimal doping point, and different properties on each 
side. For example: the underdoped samples have a lin- 
ear temperature-dependent resistivity, while the over- 
doped samples seem to act like normal metals with In- 
dependent resistivity. A pseudo-gap is found in the un- 
derdoped side, but is absent from the overdoped side. 
The magnetic penetration depth A has a minimum at op- 
timal dopingi*2i2i4. In contrast Hall measurements^ done 
in single crystals of L^-xSr^CuCU (La214), for example, 
show a smooth increase in the hole density with increas- 
ing doping. This peculiar behavior led to the general be- 
lief that the underdoped and overdoped regimes should 
be treated by different theories. While the overdoped 
side is believed to adhere to mean field theory similar to 
the BCS, the presence of strong electronic correlation in 
the underdoped side motivated the development of ex- 
otic theories for this regime. Yet, few authors found in 
the past that contrary to common belief, the hole den- 
sity is saturated in the overdoped side; for example in 
La214 6 , Bi 2 Sr 2 CaCu 2 8+(> - (Bi2212) and Tl 2 Ba 2 Cu0 6+5 
(T12201)i, and HaBa 2 Cu0 4+ 5 (Hal24)£. This possibility 
could change the way we think about overdoped cuprates. 
However, up until now, it could not be tested for the very 
popular and most homogenous YBa2Cu30 y (Y123) sys- 
tem since it could not be overdoped without introducing 
disorder, e.g., by doping with Ca. 

In this letter we overcome this problem by in- 
vestigating the fully enriched cuprate system 
(Cao.iLao.9)(Bai.66Lao.3B)OuaO y (CLBLCO). This 
system is unique in the sense that doping can vary across 
the full range, from the very underdoped to the extreme 
overdoped, without any structural changes^. The system 
belongs to the Y123 family, but it is tetragonal over the 
entire doping range, with no preferred direction for the 
CuO chains and no chain ordering. Doping is controlled 
only by the oxygen level^. We determine the doping 
level by extracting the nuclear quadrupole resonance 
parameter vq from 63 Cu Nuclear Magnetic Resonance 
(NMR) measurements. The Cu, with its spin 3/2 nuclei, 
is directly coupled to charge degrees of freedom via the 
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FIG. 1: (Color online) Symbols: The spectra of the optimally 
doped CLBLCO sample (y = 7.135). Dashed line: Powder 
line shape reconstructed using the parameters found by solv- 
ing Eqs. [I] Solid line: Best fit to the data by taking into 
account a distribution of vq. 

electric field gradient (EFG), and vq is a measure of 
this coupling, vq, in turn, depends linearly on the hole 
densitjiii n according to 

v Q = An + v Q (1) 

where A and Vq are doping- independent. This linear de- 
pendence was demonstrated for various compounds such 
as Y123 12 , La214 13 , and Hal24 14 . Therefore, Eq. [T]and 
the ability of NMR to detect the in-plane cooper [Cu(2)] 
vq selectively, will allow us to determine the carrier con- 
centration in the overdoped regime. 

The measurements were done on powder samples since 
CLBLCO is not orientable. Their preparation is de- 
scribed in Ref.-. The oxygen content was measured by 
double Iodometric titration. The accuracy of this method 
in enriched CLBLCO is about 0.01. We measured seven 
different samples in the normal state at 100 K. The most 
overdoped sample is not superconducting. The NMR 
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measurements were done by sweeping the field in a con- 
stant applied frequency / app = 77.95 MHz using a ir/2 - 
it echo sequence. The echo signal was averaged 100, 000 
times and its area evaluated as a function of field. The 
full spectrum of the optimally doped sample (y = 7.135) 
is shown in Fig. ^ Four peaks associated with the plane 
Cu(2) and one peak from the Cu(l) can be clearly seen. 

A zoom on the main features of the Cu(2) signal of 
all seven samples is depicted in Fig. |3 (note the three 
axis breakers). The evolution of the main peaks with 
doping from underdoped to the optimally doped (y — 
7.135) samples is highlighted by the dotted line. This 
line is extended to the overdoped samples to demonstrate 
their opposite peak position evolution with doping. For 
overdoped samples, arrows indicate the peak positions. 
The opposite peak position evolution of the overdoped 
samples is the main analysis independent finding of this 
work. It shows that charge density in the CuC>2 plane is 
not increasing in the overdoped side. 

The Cu spin Hamiltonian can be written aai^: 

H/h = -nl • (1 + K) ■ H+^[3I^I 2 + ?7 (I2 (2) 

where vi — ( 63 j/2tt)H, H is a unit vector in the di- 
rection of the field, K is the shift tensor, and rj is the 
asymmetry parameter of the EFG. In the absence of mag- 
netic field, there is only one transition frequency given by 
/ = VQyJ\ + ?y/3, so vq cannot be separated from 77, and 
the use of the magnetic field is essential. This field, ap- 
plied in the direction 9 and ip with respect to the principal 
axis of the EFG, lifts this degeneracy and three transi- 
tion frequencies j/ m (if, 9, ip) are expected: a center line 
which correspond to the 1/2 — > —1/2 transition (m = 0), 
and two satellites which corresponds to the 3/2 — > 1/2 
(to = 1) and —1/2 — > —3/2 (to = —1) transitions. Ex- 
pressions for v m (H, 9, ip) up to second order perturbation 
theory in vq for completely asymmetric EFG and shift 
tensors are given in Ref^&. 

For calculating a field-swept powder spectrum P(H) 
one has to evaluate 

P rn (H) ocJ2[ m2s IUp ~ "m(H, 9, tp)]dn (3) 

where M is a matrix element. This spectrum is known 
to have four peaks (two emerge from m = 0) at fields 
giveni&ii by 
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FIG. 2: (Color online) Spectra for the 7 CLBLCO samples. 
Position of the four peaks associated with the plane Cu as 
function of doping from under to optimal (y = 7.135) dop- 
ing are shown by the dotted lines. These lines are extended 
further to show the opposite peak position evolution for over- 
doped samples. Their peak position are shown by arrows. 



where vi s , vi c , Vhc, ^hs are the low satellite, low center, 
high center, and high satellite peaks, respectively, and 
we assume K x = K y — K± since CLBLCO is tetragonal. 
Usually one can extract all the averaged Hamiltonian pa- 
rameters from the position of the peaks in the spectra by 
solving Eqs.0]numericallyii. However, we are interested 
in both the parameters and their distribution. Therefore, 
we use a grid in the (9, </>) space, calculate the frequency 
v m {H, 9, ip) for every m, field, and point on the grid, and 
add 1 to a histogram of H when this frequency equals 
/app- The matrix elements are taken as unity. The nu- 
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merical summation approximates P(H) in Eq. [3] 

This numerical approach is demonstrated in Fig. ^ 
where the four peaks associated with Cu(2) have almost 
axial symmetry (rj ~ 0), and the fifth (middle) peak 
is from Cu(l) which has a completely asymmetric EFG 
(r) ~ 1). The dashed line in Fig. ^ is the powder line 
shape calculated after the parameters for Cu(2) were ex- 
tracted by solving Eqs. 0] for the plane site. The sin- 
gularities of the theoretical line agree with the peaks in 
the NMR data, as expected. However, the overall shape 
of the calculated line does not resemble the data at all, 
the reason being the distribution of the Hamiltonian pa- 
rameters. The main contribution to the width, of the 
order of a few MHz, is from a distribution in vq, since 
the quadrupole interaction is the only interaction of such 
magnitude in the system. To account for this distribu- 
tion we simulated the line shape by adding to the nu- 
merical evaluation of Eq. [3] a loop over 200 values of vq 
drawn from a normal distribution with a width Avq . We 
also added the contribution of the chain site, with rj ~ 1 
and removed the K x = K y = K± constrain, to take 
into account possible local orthorhombic distortions. We 
searched for the best fit to the data by x 2 minimiza- 
tion using a simplex codsi&. The result for the optimally 
doped sample is shown as the solid line in Fig. ^ 

In Fig.[3]we show the position of the four peaks (shown 
in Fig. ^ as a function of doping. The separation be- 
tween the satellites increases with doping up to the op- 
timal doping point and decreases as the system is over- 
doped. The evolution of the width of the lines can be 
seen by examining the peaks of the center lines in Fig El 

By solving Eqs.^J we calculate i/q, t], K z and K± for 
all the samples. The calculated vqs are shown in Fig. 

All parameters are used as the initial guess in the fit- 
ting program. The results for fitted vq and Avq are also 
shown in Fig. 0] The solid lines in this figure are guides 
for the eye. The width of the lines does not allow us to 
determine rj, K z , K x and K y very accurately. In gen- 
eral the Cu(2) rj is always very small for all the samples, 
ranging from 0.015 for the optimal doped sample and up 
to 0.07 for the underdoped sample, in agreement with 
Y123i£. We have also attempted to fit the data starting 
from arbitrary initial guess. This gave somewhat differ- 
ent fit values. The error bars in Fig.QJare estimated from 
the different fit procedures. These error bars are larger 
by a factor of 4 than the error bars calculated by the 
fitting program itself. 

The location of the simulated peaks agrees with the 
data for all samples, but the simulation program yields a 
higher vq than extracted from Eqs.0] The reason is that 
the theoretical peaks for zero Avq are not symmetric (see 
dashed line in Fig. QJ. By introducing a normal distri- 
bution of vq the peaks tend to shift towards the center 
of the line resulting in what seems to be a smaller aver- 
age vq. The computer fit program corrects for this effect 
by selecting a higher average vq than the initial guess 
(based on Eqs.QJ. We believe that the fit results repre- 
sent the data better than the calculation that is based on 
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FIG. 3: (Color online) Position of the 4 peaks associated with 
the plane Cu as function of doping. 



the peaks alone. 

The results indicate that vq grows linearly with doping 
in the underdoped side of the phase diagram in agree- 
ment with other compoundsiSii&H. In the overdoped 
side vq grows very slowly with doping and saturates; it 
does not decrease as one might think looking at the peaks 
position only. This is due to the increase in Avq. The 
width decreases with doping, having a minimum at opti- 
mal doping, and increases quite sharply in the overdoped 
side. We emphasize, that the main result, a deviation 
from the linear dependence of vq in the overdoped side, 
is observable in the raw data itself. The same result is 
supported both by the exact solution of Eqs. 0] and by 
the numerical fits. 

The fact that there is a sudden change in the vq ver- 
sus y relation suggests a change in the doping mechanism 
itself. In CLBLCO, like in Y123, the doping is done by 
adding oxygen ions into the plane of chains, which are 
initially only half full. Our result suggests that in the un- 
derdoped region a constant part of the introduced holes 
goes into the Cu02 plane. After crossing the optimal 
doping point, most of the holes remain in the chains. 

Recent density functional calculations of doping and 
its effect on the EFGSi in YBCO show that both vq and 
the doping should go through a maximum as a function 
of the oxygen level. It should be noted that resistivity 
and Hall measurements in YBCO will be less sensitive to 
the change in the plane doping mechanism since above 
optimal doping the chain layers start to conduct. 

Thus, above optimal doping the energetic cost of 
adding a free hole to the Cu02 plane increases. Trying 
to overdope CLBLCO results in transfer of charge into 
the chains. Moreover, attempts to overdope the sample, 
leads to higher degree of disorder in the plane, manifested 
in the increase of Avq. This is the main conclusion of 
this work. 
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FIG. 4: (Color online) (a) vq, the mean value of the quadru- 
ple frequency, versus doping as obtained by fit and calculation 
described in the text, (b) Avq, the width of the normal dis- 
tribution around vq obtained by the fit. 

Finally, it is interesting to compare the NQR re- 
sults with previous magnetic penetration depth, A, 
measurements^. A was determined by transverse field 
muon spin rotation where the relaxation rate a is pro- 
portional to A~ 2 [ref 21 ]. In CLBLCO a is proportional 



to T c in both the under and over doped sides, with the 
same proportionality constant of all HTSC, obeying the 
Uemura relation^. Within the clean limit BCS theory 
framework, at T — all the normal state carriers should 
pair and condensate giving tr oc n s /m* where n s = n/2 
and the effective mass is just twice the effective mass of 
the quasi-particle in the normal state. Therefore, the sat- 
uration of the number of holes in the plane above optimal 
doping cannot explain the reduction in a alone. Clearly 
the explanation must involve a combination of the satu- 
ration of the normal state carrier density and some other 
mechanism. Few suggestion for the nature of this mech- 
anism were proposed, for example, phase separation and 
creation of metallic droplets^, and pair breaking 2 . The 
fact that we find an increase in the width of the charge 
distribution in the planes might be an important ingredi- 
ent in the correct theory explaining the peculiar behavior 
of the superfluid density in the overdoped side. 

In summary, we have managed to measure the average 
quadrupole resonance parameter vq and the width of its 
distribution Avq, for the in-plane Cu, in a system be- 
longing to the YBCO family but where doping can vary 
from heavily underdoped to extreme overdoped. We find 
that vq increases with oxygen doping in the underdoped 
side but is saturated in the overdoped side. In contrast, 
Avq has a minimum at optimal doping. 

We acknowledge useful discussion with Y. Eckstein and 
the help of Arkady Knizhnik with the Iodometric titation. 
This work was funded by the Israeli Science Foundation. 



* Current Address: Department of Physics, University of 
Illinois at Chicago, IL 60607 

1 Y.J. Uemura, A. Keren, L.P. Le, CM. Luke, W.D. Wu, Y. 
Kubo, T.Manako, Y. Shimakawa, M. Subramanian, J. L. 
Cobb and J. T. Markert Nature, 364 605 (1993). 

2 Ch. Niedermayer, C. Bernhard, U. Binninger, and H. Glck- 
ler Phys. Rev. Lett., 71 1764 (1993). 

3 J. P. Locquet, Y. Jaccard, A. Cretton, E. J. Williams, F. 
Arrouy, E. Mchler, T. Schneider, O. Fischer and P. Marti- 
noli Phys. Rev B, 54, 7481 (1996). 

4 A. Keren, A. Kanigel, J. S. Lord and A. Amato Solid State 
Comm., 126 39 (2003). 

5 Yoichi Ando, A. N. Lavrov, Seiki Komiya, Kouji Segawa, 
and X. F. Sun Phys. Rev. Lett., 87, 017001 (2001). 

6 S. Uchida, T. Ido, H. Takagi, T. Arima, Y. Tokura, and S. 
Tajima Phys. Rev. B, 43 7942 (1991). 

7 A. V. Puchkov, P. Fournier, T. Timusk, and N. N. 
Kolesnikov Phys. Rev. Lett, 77, 1853 (1996). 

8 R. Puzniak, R. Usami and H. Yamauchi Phys. Rev. B 53, 
86 (1996). 

9 D. Goldschmidt, A. Knizhnik, Y. Direktovitch, G. M. Reis- 
ner, and Y. Eckstein Phys Rev. B 49, 15928 (1994). 

Y.J. Uemura Solid State Comm. 120 , 347 (2001). 

1 Kunisuke Asayama, Yoshio Kitaoka, Zheng Guo-qing and 
Kenji Ishida Prog. Nuc. Mag. Reso. Spec, 28, 221 (1996). 

2 HYasuoka, in Spectroscopy of Mott Insulator and Corre- 



lated Metals, edited by A. Fujimori and Y. Tokura, Solid 
State Sciences, Vol. 119 (Springer- Verlag, Berlin, 1995), p. 
213. 

13 Guo-qing Zheng, Yoshio Kitaoka, Kenji Ishida and Ku- 
nisuke Asayama, J. Phys. Soc. Jpn. 64, 2524 (1995). 

14 A. A. Gippius, E.V. Antipov, W. Hoffmann and K. Luders, 
Physica C 276, 57 (1997). 

15 CP. Slichter Principles of Magnetic Resonance, Harper 
and Row, New York, (1963). 

16 J. F. Baugher, P.C Taylor, T. Oja and P.J. Bray J. of 
Chem. Phys., 50, 4914 (1969). 

17 A. Keren, P. Mendels, M. Horvati?, F. Ferrer, Y. J. Ue- 
mura, M. Mekata and T. Asano Phys. Rev. B 57, 10745 
(1998). 

18 W.H. Press, B.P. Flannery, S. Teukolsky and W.T. Vetter- 
ling, numerical Recipes in C, Cambride University Press 
(1988). 

19 C. H. Pennington, D. J. Durand, D. B. Zax, C. P. Slichter, 
J. P. Rice, and D. M. Ginsberg, Phys. Rev. B 37, R7944 
(1988); T. Shimizu et al, Journal of the Physical Society 
of Japan, 57, 2494 (1988). 

20 E. P. Stoll and P. F. Meier, Private Communication. 

21 Muon Science: Muons in Physics, Chemistry and Materi- 
als, Eds S. L. Lee, S. H. Kilcoyne and R. Cywinski (Insti- 
tute of Physics, London), 1999. 



